Mott insulators in plaquettes 
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We study small systems of strongly interacting ultracold atoms under the influence of gauge 
potentials and spin-orbit couplings. We use second order perturbation theory in tunneling, derive 
an effective theory for the strongly correlated insulating states with one atom per site, and solve 
it exactly. We find dramatic changes in the level structure and in the amount of degeneracies 
expected. We also demonstrate the dynamical behavior as the barriers between plaquettes are 
gradually removed and find potentially high overlap with the resonating valence bond (RVB) state 
of the larger system. 
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I. INTRODUCTION 

Weakly interacting gas of bosons can form a Bose- 
Einstein condensate at low enough temperatures. This 
condensate can be often well described by a classical the- 
ory where condensate atoms are described by a conden- 
sate wavefuction whose dynamics follows from the Gross- 
Pitaevskii equation [lj. I n optical lattices potential bar- 
riers between sites suppress the movement of atoms and 
increase confinement which increases effective interaction 
strengths. Consequently optical lattices can push the 
weakly interacting quantum system into the regime a 
strongly correlated physics. Bose-Einstein condensation 
of weakly interacting systems gives way to Mott-insulator 
physics at high interactions and this transition was ex- 
perimentally observed by Greiner et al. @- 

While the position of the transition (in higher dimen- 
sional systems) can be fairly accurately predicted with 
the Gutzwiller ansatz that is a product state of individual 
site quantum states, it predicts the physics of Mott insu- 
lators poorly since it does not take into account higher 
order tunneling processes. However these strongly corre- 
lated insulators can be studied in terms of different spin 
models and can have a very rich physics of quantum mag- 
netism and criticality Q. At a formal level spin models 
can be derived from the second order perturbation theory 
in kinetic (tunneling) energy of the lattice model 0, Q . 

Recently there has been rapid experimental progress in 
creating artificial spin-orbit couplings @ as well as gauge 
potentials 0, HI in ultra cold quantum gases. Given this 
background even non-Abelian gauge fields [1, have 
become feasible. _[For a recent review see, for example, 



Dalibard et al. In addition to these advances, a 



three-dimensional array of independent plaquettes has 
also been demonstrated even to a level of being able to 
prepare specific quantum states in each plaquette 12]. 

The purpose of this article is to study insulating states 
in an optical lattice for small systems by exactly diago- 
nalizing the effective Hamiltonian corresponding to the 



second order perturbation theory in tunneling. We will 
solve the system under quite general conditions so that 
tunneling coefficients can be spatially varying and arti- 
ficial gauge potentials included easily. We demonstrate 
the level structure for few simple choices of phases and 
point out when, for example, the amount of degeneracies 
can vary with gauge potentials. Furthermore, we show 
some examples of dynamics in multi-plaquette systems 
and solve the problem of plaquettes with spin-orbit in- 
teractions exactly for 2- and 4-sites. 



There are some interesting recent studies on the phase 
diagrams of Mott insulators with spin-orbit couplings [HI, 
[HI that pointed out how the effective Hamiltonian be- 
comes a combination of Heisenberg model, quantum com- 
pass model, and Dzyaloshinskii-Moriya interactions. The 
approach used in this paper is different in that our focus is 
on the full solution for small systems which naturally en- 
able also the solution of time-dependent problems. This 
is relevant since in the experiment by Nascimbene et 
al. [l2j resonating valence bond states were prepared in 
plaquettes of an optical lattice. In these experiments dy- 
namics was an important ingredient and since the system 
was made out of independent plaquettes, exact diagonal- 
ization for small systems provides a natural theoretical 
framework. 



' Jani-Pctri.Martikaincn@aalto.fi 



The paper is organized as follows. We begin by out- 
lining the theory relevant for our purposes in Sec. |IIJ In 
Sec. lIIII we exactly diagonalize the system with some sim- 
ple gauge potentials both for the square as well as for the 
triangular plaquette. Then in Sec.|lV]we explore the dy- 
namical behavior when barriers between plaquettes are 
removed and in Sec. |V]we exactly diagonalize the Mott 
insulating plaquette with spin-orbit interactions. We end 
with a few concluding remarks in Sec. I VII 



II. FORMALISM 



In each site r of the lattice the local physics for bosons 
is described by the interaction Hamiltonian [H, [f| 



#J,r = -^-^l,r(«l,r - 1) 

+ U^flXrfll.r, 



(1) 



where n a , T = ?/>J, r ip a ,r (a <E {1, 2}) are the density opera- 
tors for each component. Here the onsite eigenstates can 
be written in terms of occupation numbers as \ni,n2) r - 
The approach works equally well for fermions, but then 
Uu = U22 — since s-wave interaction between identical 
fermions vanishes and Pauli exclusion principle must be 
enforced. Furthermore, this problem can be straightfor- 
wardly generalized to spinorial bosonic systems, multi- 
component fermionic systems [15| . or Bose- Fermi mix- 
tures [HI [TtJ by s imp ly changing the term describing 
onsite interactions (l8| |. 

In this work tunneling processes are kept very general 
so that they can include position dependent tunnelings 
as well as, later in the text, spin-orbit couplings. In the 
absense of spin-orbit couplings the tunneling processes 
are described by 



1 V^a, 



(2) 



where a indicates the component index, m indicates the 
bond, while ipa,m+ and ipa,m- are the field operators for 
atoms of type a on either end of the bond m. Later we 
will introduce also spin-orbit couplings in which case the 
atom type can change in the tunneling process. 

Throughout this paper we focus on the experimentally 
simplest case of having one atom per site. For all states 
having just one atom per site interaction energy vanishes 
so double (and higher) occupancy is energetically unfa- 
vorable. Adding tunneling terms into the model implies 
removing an atom from one site and adding it into an- 
other. If the initial state only had singly occupied sites, 
this implies that tunneling process creates a state where 
one site is doubly occupied. In the lowest order pertur- 
bation theory the tunneling term does not give rise to 
any energy shift since the state with double occupancy 
is orthogonal to initial state, but at second order there 
can be another tunneling process which brings the state 
back into the initial subspace of singly occupied sites. At 
this level tunneling terms do matter and they determine 
the energetics and dynamics of the various possible Mott 
insulating states. 

When the tunneling processes are described by a 
Hamiltonian Ht this second order perturbation theory 
gives rise to an effective Hamiltonian [lil ] 



E 7 -(E a + E )/2' 



(3) 



where E{ are the eigenenergies of the interaction Hamil- 
tonian, a and (3 live in the manifold with one atom per 
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FIG. 1. (a) Schematic description of the 2x2 plaquette and 
notation used in the text, (b) Schematic description a trian- 
gular plaquette. 



site, while states labeled by 7 have one site with 2 atoms 
per site. 

In the literature the this expression is typically writ- 
ten in terms of pseudo-spin operators. While this can 
formally be done in our cases as well, the expressions be- 
come ever more complex as tunneling properties in the 
system become more diverse. For this reason we choose 
not to rephrase the problem in terms of spin-operators, 
but construct the relevant Hamiltonians in the subspace 
with one atom per site with the computer. In this way 
all associated "spin" -models are included, but as an in- 
put one only needs the onsite interactions and tunnel- 
ing coefficients. In principle, both interactions as well as 
tunnelings could be position dependent. The latter op- 
tion is in fact important when we wish to study effects 
of gauge potentials in lattices and position dependence 
of atom-atom interactions can be important, for example, 
in bipartite lattices fl9j . If we were to write results in 
terms of pseudo-spins, we would have to derive a sepa- 
rate spin-model for each variant of the problem. 

Our procedure amounts to first constructing all eigen- 
states of H[ = J2 r Hi,* w i tn ( on average) unit filling 
and at maximum 2 atoms per site. From these states we 
construct the subspace with exactly one atom per site 
and then compute the elements H e ff numerically. The 
resulting Hamiltonian can then be diagonalized and this 
solution (if so desired) can then be used to solve dynam- 
ical problems as well. In this paper we will focus on 
2 x 2 or triangular plaquettes since the first system has 
been experimentally realized [l2[ and the latter one can 
be realized [20| . Furthermore, both can be solved fully 
without difficulty. We also solve dynamical problems in 
systems of two and three weakly coupled plaquettes. 



III. MOTT INSULATORS IN GAUGE 
POTENTIALS 

We show a schematic description of our system(s) ele- 
mentary cells in Fig.[TJ The various tunneling coefficients 
are given by t m and they describe processes like (exam- 
ple for ti) -t a ^ a 2 ^ a ,i - t* a .i4>i A 4> a ,2- for tunneling of 
a atoms between sites 1 and 2. Nascimbene et al. fijl 
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solved the eigenstates of this system in the case when 
t 2 = ti = t y is different from t\ = t 3 = t x . We wish 
to understand the role of gauge potentials in such a sys- 
tem and for this reason we allow for the possibility of 
complex and position dependent tunneling coefficients. 
In contrast to superfluid regime, if both components ex- 
perience the same gauge potential, gauge potential does 
not affect the Mott insulator. For gauge potentials to 
play a role, we have to allow for the possibility of spin 
dependent gauge potential. To keep things as simple as 
possible we choose tx,i — all other tunneling co- 

efficients as 1/2, and vary the phase 0. (Using Peirls 
substitution the phase could be related to gauge vector 
potential A through (f> — J r 3 dr ■ A/h). This gives rise 
to the Harper model in a plaquette [2l| and if the flux 
a — (p/2ir = p/q is rational the ideal theory has a g-fold 
degenerate ground state. 

In Fig. [2] (a) we show the eigenenergies as a function 
of a. At a = the result agrees with those reported 
in Ref. [l2j], but as gauge potential is turned on the 
3-fold degeneracy of the fist excited state is broken as 
the d-wave RVB state splits into a separate branch. As 
gauge potential becomes stronger the ground state be- 
comes doubly degenerate at a = 1 /2 also in the interact- 
ing system. Furthermore, the gap separating the highest 
s-wave RVB state closes at a = 1/2 and the eigenstate 
there becomes doubly degenerate. 

One way to create gauge potentials in lattices is to 
shake the lattice [13, appropriately and such tools 
have indeed been demonstrated experimentally [H, l23l - 
[25j . Some of the more simple methods of shaking can- 
not create a magnetic flux in a square lattice while they 
might do so in a triangular lattice. For this reason it is 
instructive to solve the problem also for a triangular pla- 
quette. We show a typical the result in Fig. [2] (b). In a 
triangular plaquette with one atom per site the different 
component have necessarily different atom numbers. If 
N% = N2 + 1 there are three eigensolutions. The excited 
state is doubly degenerate in the absence of gauge poten- 
tial while its degeneracy is broken in the gauge potential. 
At a — 1/2 the ground state is again doubly degenerate 
while the excited state is non-degenerate. 



IV. FEW COUPLED PLAQUETTES 

In the experiment by Nascimbene et al. [l2T | the pla- 
quettes were not strongly coupled and they were able 
to create the s-wave RVB state in each plaquette. We 
can use our approach to also study systems with more 
than one plaquette. In this case it becomes possible to 
study how coupling between plaquette changes the sys- 
tem behavior. Since we have not made any assumptions 
about the magnitudes of tunneling strengths we can sim- 
ply add tunneling terms between plaquettes without, for 
example, a need for a new perturbation theory in "small" 
inter-plaquette tunneling. One notable result from solv- 
ing the problem for coupled plaquettes (for bosons) is 
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FIG. 2. (a) Eigenenergies of 2 x 2 plaquette when another 
component experiences a gauge potential and iVi = N-2- The 
uppermost state at a = is the s-wave RVB state and the 
first excited state is the d-wave RVB state as discussed in 
Ref. [H. We choose \U\ = 1/2 and Un = U22 = U12 = 
1. (b) Eigenenergies in a triangular plaquette when another 
component experiences a gauge potential and Ni — N2 + 1. 



that the ground state is degenerate in the limit of van- 
ishing inter-plaquette tunneling and is separated from the 
first excited states by only a small gap as inter-plaquette 
tunneling is introduced. 

In sharp contrast to this, the most excited state cor- 
responding to the product state of plaquette s-wave 
RVB states is gapped in the limit of vanishing inter- 
plaquette tunneling and the gap remains quite large 
even as the inter-plaquette tunneling approaches the in- 
plaquette tunneling. More quantitatively, when all cou- 
pling strengths are the same, Un = U22 = U12 = U, 
and in-plaquette tunneling is taken as t = 1/2 the gap 
AE = l/U when inter-plaquette tunneling is zero. When 
inter-plaquette tunneling equals in-plaquette tunneling 
the gap is reduced to AE ~ 0.77/J7 for two plaquettes 
and AE 0.65/C/ for three. This separation of the (ex- 
perimentally relevant) most excited state might simplify 
studies where Hamiltonian is changed as a function of 
time, since unwanted transitions to other states can be 
suppressed. 

To understand coupled plaquettes better, let us then 
demonstrate the dynamical behavior for two coupled pla- 
quettes. We take the initial state to be the most excited 
state for uncoupled plaquettes. This state corresponds 
to the product state of plaquette s-wave RVB states. We 
then either increase the inter-plaquette tunneling tb by 
linearly ramping it to the same value as the in-plaquette 
tunneling t or by turning it on instantaneously. In Fig. [3] 
we show how overlaps with the initial state and the RVB- 
state over all 8 sites (equal superposition of all possible 
ways to cover the lattice with singlet bonds between near- 
est neighbors) behave. As can be seen, the initial state 
has a substantial overlap with the resonating valence 
bond state in the 8-site system and this overlap increases 
further as the inter-plaquette tunneling approaches the 
in-plaquette tunneling and bonds start forming also be- 
tween plaquettes. In this example, the maximum over- 
lap is around 97%. The overlap with the initial prod- 
uct state remains large, but is substantially smaller than 
the overlap with the RVB-state. Incidentally the overlap 
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FIG. 3. Overlap O = \{RVB\ip(t))\ 2 with the RVB-state over 
all 8-sites (solid line) and with the initial state \ipo) that was 
a product state of individual plaquette s-wave RVB-states 
(dashed line). Insets shows the strength of the inter-plaquette 
tunneling tb relative to in-plaquette tunneling t. In (a) we 
choose to ramp up the inter-plaquette tunneling strength lin- 
early and the kept it steady for the remainder of the sim- 
ulation. In (b) the inter-plaquette tunneling was abruptly 
turned on to be equal to in-plaquette tunneling. We used 
Uii = U22 = C/12 = 1 and t — 1/2, while the unit of time was 
t = 1/U 11 . 



of the RVB-state with the most excited state of the 8- 
site system when inter-plaquette tunneling is equal to in- 
plaquette tunneling is around 96%. Our result with the 
linear ramp demonstrates that the most excited state of 
the uncoupled plaquettes evolves adiabatically very close 
to the most excited state of two strongly coupled plaque- 
ttes. 

In Fig. |3] (b) we show that if the barriers between pla- 
quettes are removed suddenly the system, unexpectedly, 
strays further away from the initial state and dynamical 
behavior is more dramatic. These results demonstrate 
that at least for small systems RVB-states can be formed 
with fairly high fidelity starting from a collection of RVB- 
states in independent plaquettes and lowering the barri- 
ers between plaquettes. 



V. PLAQUETTES WITH SPIN-ORBIT 
COUPLING 



Recently it has become possible to create spin-orbit 
couplings in ultracold atomic gases @. This exciting de- 
velopment has opened up entirely new possibilities and 
motivated a rapidly expanding literature on various as- 
pects of spin-orbit coupled systems using cold gases. The 
presence of spin-orbit coupling (SOC) changes our pre- 
vious discussion in that it becomes possible for atoms to 
change type as they tunnel between sites. The tight bind- 
ing model for a spin-orbit coupled two-component gas is 
given by € {1,2}) 



EE 

rij ~f—x,y 



+ H.c. 



(4) 



where ip. 
and t? 7 



t 



creates an atom of type i in the lattice site r 
07 is the vector to the nearest neighbor of r 



along the direction 7 (a is lattice spacing). T^lj^ are the 
tunneling matrices, but note that (as in previous section) 
since tunneling coefficients in the plaquette can be made 
to vary, we have kept the possibility of position depen- 
dence in them. In our notation we mostly follow Radic 
et al. [3]. 

The free spin-orbit coupled system can be expressed in 
terms of a vector potential A with a Hamiltonian Hq = 
(p — A) 2 /2m, where m is the atomic mass. We choose 
the vector potential to be position independent and given 
by 



A = (— maa x , —m/3a y ) , 



(5) 



where <r 7 (7 € {x, y, z}) are the Pauli matrices. Peirls 
substitution 



—iaA-, 



(6) 



can then be used to relate tunneling coefficient t~ in the 
absence of SOC (i.e. 9 X = 9 y = 0) to those with SOC. 
This substitution is not exact and can be improved, but it 
is reasonably accurate over wide range of parameters and 
where its accuracy is somewhat worse it is still qualita- 
tively useful (l4| . Inaccuracies can become more serious 
in and closer to the superfluid regime. 

When the expressions for tunneling elements in SOC 
system are expanded we find 



j, _ , ( cos 9 X i sin X 
x ~ x I i sin 6 X cos 6 X 



and 



1 7, 



cos 6 y sin 6 y 



(7) 



(8) 



In Fig. [4] we show the eigenenergies of the spin-orbit 
coupled system with only 2-sites. The basis is spanned 
by 4 states and without spin-orbit coupling the ground 
state in 3-fold degenerate. Spin-orbit coupling breaks 
this degeneracy so that ground state remains 2-fold de- 
generate. These ground states do not change with spin- 
orbit coupling, but the 2 excited state energies do de- 
pend in the strength of SO coupling. The eigenstates 
are however always such that the lower one is of type 
|V>) = (|10)i|10) 2 -|01)i|01) 2 )/V2 while the upper one is 
of type |V) =(|10)i|01) 2 -|01)i|10} a )/V2. 

In a two-dimensional system situation is more complex. 
In Fig. [S] we show the eigenenergies of the spin-orbit cou- 
pled plaquette as a function of 9 X and 9 y . As SO coupling 
is turned on the initial three-fold degeneracy of the first 
excited state is broken and all states are non-degenerate. 
However, remarkably as one approaches 6 X = 9 y = n/2 
all states become degenerate. (For non-interacting sys- 
tem one has at this point three-fold degenerate states 

with energies E± = ±^Jt x + t y .) This suggests that close 

to this point the perturbative approach we have used 
might be breaking down and more accurate theory might 
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FIG. 5. (Color online) Eigenenergies for a 2 x 2 system with 
spin-orbit coupling as a function of 9 X and 9 y . The inset in 
(a) shows the solutions at 9 y = 0. We choose Uu = U22 = 
U12 — 1 and t x = t y = 0.5. All states are degenerate at 
9 x = 9 y = tt/2. 

be required to properly resolve the eigenstates. Also, this 
high dimensionality of the ground state manifold might 
have interesting consequences for the dynamics as well as 
on how the system responds to experimental probes. 
VI. CONCLUSIONS 

In this article we have explored the physics of small 
insulating plaquettes under the influence of gauge po- 
tentials and spin-orbit couplings. Computations were 
done for bosonic systems, but in the simplest cases two- 
component fermionic problem amounts to sign change in 
the eigenenergies. For example, the highest excited state 
for bosons can become the ground state for fermions. At 
the formal level our approach does not distinguish be- 
tween bosons and fermions since the only inputs into the 
solver are the various onsite interactions strengths and 
tunneling coefficients. For this reason the approach used 
here is also easy to apply to the study of fermionic multi- 
flavor systems [l5| or Bose-Fermi mixtures. Since atoms 
could also be prepared on the excited bands of the lat- 
tice [26{ , it would be of interest to also explore the Mott 
insulating plaquettes in bipartite systems |19l | and in ex- 
cited bands. 
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